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I. INTRODUCTION 



Thin films of nematic liquid crystals (NLC) have attracted attention over the years, as evidenced 
by a number of experimental and theoretical studies^^°. When thin nematic films are deposited 
on solid or liquid substrates, they often exhibit antagonistic anchoring at the free surface and at 
the substrate, i.e., the director orientation at the substrate is generally parallel to the substrate 
(planar anchoring) but at the free surface the director is orthogonal to the surface (homeotropic 
anchoring). As a consequence, the local director orientation changes across the film resulting in an 
elastic contribution to the energy that should not be neglected: such films are called hybrid films. 
Sometimes instabilities are observed that result in lateral periodic stripe patterns of the director 
orientation'^ and film height. However, this is only the case for thin films with thicknesses of 
several hundred nanometer and below; the wavelength of the stripe patterns diverges at an upper 
critical film thickness and so, for thicker films, only the usual defects of the nematic phases are 
observed^'^. Note that spinodal patterns have also been observed**'^' normally, in the vicinity of 
nematic-isotropic or smectic-nematic phase transitions. In contrast, the stripe patterns are observed 
well inside the nematic region of the liquid crystal phase diagram. 

In order to develop a theory for the behaviour of confined nematic liquid crystals, one may 
calculate the director orientation profile for a given static free surface. Typically, either a flat 
film or a periodically deformed state is considered. Such a given static geometry is then used to 
investigate the director field and to determine its stability. For an imposed flat film, an energy 
argument allows one to show that there exists a critical thickness'"* 

K K 

he = r- 

A+ A_ 

where K is the bulk elastic constant of the liquid crystal (in the one constant approximation) 
and and A_ are the anchoring strengths at the free surface and at the substrate, respectively. 
For thin films, with thickness h < he, the director profile is undistorted; the film is in the so- 
called planar (P) state and the director is aligned parallel to the anchoring angle at the interface 
with the stronger anchoring strength. For thick films, with film thickness h > the state that 
minimises the free energy is that where the director orientation changes continuously between 
the two anchoring directions as one moves across the film; this is the Hybrid-Aligned-Nematic 
(HAN) state introduced above and is the case for the strong anchoring situation considered here. 
If one assumes that the system is invariant in one direction across the surface on which the film 
is deposited, so that the film is effectively two-dimensional (2D), one finds that these states are 



(1) 



linearly stable. To confirm this assumption, much effort has gone into determining whether the 
film is laterally stable'"^''°. However, since the film geometry is imposed and static, such analyses 
can not account for a possible coupling of variations in film height and director orientation. 

In alternative approaches, film thickness evolution equations for films of nematic liquid crystals 
are derived employing lubrication or long-wave hydrodynamic theory. This approach is used 
successfully to explore the dynamics of thin films of simple liquids under the influence of gravity 
or other body forces, and a variety of surface and interfacial forces'^ Based on the long-wave 
approach, Ben Amar and Cummings^ derived a model to describe the surface evolution of NLCs 
with strong anchoring in 2D settings that was later adapted to model 2D spreading droplets^^, 
spreading droplets with defects^"* and to account for three dimensional settings (SD)^"*. Another 
long-wave model was introduced by Carou et al. to study blade coating and cavity filling flows 
of NLC in 2D^^"^^. However, neither of these long- wave evolution equations agrees with models 
that use energy arguments, when it comes to ascertaining the effect of the elastic distortion energy 
on the film dynamics. Antagonistic anchoring is predicted to destabilise the film in Refs. 7, 22- 
24, but in Refs. 25-27 it is predicted to have no influence on the stability of the film. In Ref. 9, 
however, it is argued on physical grounds that the elastic energy is stabilising. 

The purpose of this note is to clarify this issue by consolidating the hydrodynamic long-wave 
and energetic approaches in the case of layers of nematic liquid crystals with strong anchoring, and 
so to provide a solid basis for further investigations. The manuscript is organised as follows: The 
continuum theory of NLC, including the elastic energy and Ericksen-Leslie equations together 
with consistent boundary conditions, is given in Section II. Focusing on the 2D case, the long- 
wave approximation of the governing equations and boundary conditions is sketched in Section III 
while the full details are given in Appendix A. This allows the reader to easily reproduce our 
main findings. In Section IV, a thermodynamically motivated gradient dynamics formulation is 
employed to derive the evolution equation of a nematic film. The stability of the free surface is 
studied through a linear stability analysis. Finally, in Section V we compare the results of the two 
approaches and discuss the validity and limitation of the present model. The note concludes with 
an outlook on related problems that could be studied based on our results. 
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II. CONTINUUM DESCRIPTION OF NEMATIC LIQUID CRYSTAL 



Nematic liquid crystals consist of rod-like molecules that have no positional order, but have 
long-range orientational order. Thus, the molecules are free to flow as a liquid, but still main- 
tain their long-range directional order. The mean molecule alignment is described by the unit 
vector n = (rii, n2, n^^Y where the superscript T denotes matrix transposition. Further notation 
conventions used here are presented in Appendix C. 

Distortions of the director field result in a contribution to the free energy, that for NLC is known 
as the Frank-Oseen elastic energy and reads^^'^^ 

+ ^{K2 + K^)V ■ ((n ■ V)n - (V ■ n)n) , (2) 

where Ki, K2 and are the splay, twist and bend elastic constants, respectively, and {K2 + K^) 
is called the saddle-splay constant. Note that the saddle-splay term is often omitted since it does 
not contribute to the governing equations for problems invoking strong anchoring. 
We use the one-constant approximation to simplify the problem. One assumes^^'^^ 

K = K, = K2 = K3, = 0. (3) 

and obtains the simplified energy density 

K K 
Wf = — Vn: {Vnf = — ni^kni,k (4) 

that enters the nemato-hydrodynamic equations discussed next. 



A. Ericksen-Leslie equation 

The flow of NLC may be described by the Ericksen-Leslie equations^^^^ The flow is incom- 
pressible, satisfying 

V ■ = 0, (5) 

where v = (f 1, f 2, fs)^ is the velocity field. The momentum balance equation is 

D 

p—v = V ■ cr, (6) 
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where p is the density, D / Dt = d / dt + v ■ V is the material derivative, t is the time variable and 
cr is the stress tensor of the NLC. The stress tensor is defined as^^ 

o- = -p/ + cr^ + cr^, (7) 

where p is the pressure, / is the identity tensor, cr^ is the elastic (Ericksen) stress tensor, defined 
by 

(t"" = -—f-i\/nf, (8) 
oVn 

and cr^ is the viscous stress tensor, written as 

[cr^] = ainkUpekpUirij + a2Ninj + asNjUi + a^Cij + a^eikUkUj + aQejknkfii, (9) 

where 

eij = ^ {vij + Vj^i) , Wij = ^ {vij - Vj^i) , Ni = —rii - WikUk. (10) 

The ftj are constant viscosities. 

The equation for the balance of angular momentum is written as (neglecting director inertia) 

\dVn J dn 

where 

[9]i = -liN-i - -f2ei,knk, 7i = "s - "2, 72 = - "s- (12) 

Furthermore, A is the Lagrange multiplier ensuring \n\ = 1. 

Under the assumption of the one constant approximation, the Ericksen-Leslie equations, 
Eqs. (5), (6), and (11) simplify to 

V ■v = 0, (13a) 
p^v = -V{p + wp) - KVn ■ An + V ■ (7^, (13b) 
KAn + g = Xn, (13c) 

respectively, where we have used that 

(T^ = -KVn- {Vnf, (14) 

and 

V ■ (Vn ■ (Vnf) = (Vn : (Vn)^) + Vn ■ An. (15) 

As a result, the Ericksen-Leslie equations in the one constant approximation are given by Eq. (13) 
and need to be solved subject to appropriate boundary conditions. 



a. Remark 1: Note that sometimes the stress tensor for NLC is written differently from 
Eq. (7), e.g., Ref. 32 uses 

a- = -(p + wf) / + cr^ + cr^. 

However, one may combine the two terms of the isotropic part of a and define a modified pressure 
as p = p + wp- Hence, with the exception of the modified pressure our derivations that follow are 
not affected. 

b. Remark 2: Equation (13b) can be rewritten as 

p^i; = — V(p + wp) + Vn ■ g + V ■ CT^ 

by using Eq. (13c) together with Vn • n = 0. This formulation is more popular in the literature 
since it only involves the first derivative of the director field. 

1. Boundary conditions 

We assume here that the NLC film sits on a solid substrate at 2; = with the free surface (or 
film thickness) described by 2; = h{x, y, t). 

For the director field n, we impose strong anchoring conditions such that the director is planar 
at the solid substrate and is homeotropic at the free surface. Specifically, we have 

n ■ z = 0, atz = 0, (16a) 
n • tj = 0, at z = h{x,y,t), (16b) 

where z = (0, 0, 1)-^ and ti are the surface tangent vectors, 

tl = ^ (1, 0, dM^, t2 = ^ (0, 1, dyhf. (17) 

For the velocity field v, we assume no-shp at the solid substrate, 

= V2 = vs = 0, atz = 0. (18) 

At the free surface, z = h{x, y, t), we have the kinematic condition and balance of normal and 
tangential stresses. The kinematic condition is 



W3 = dth + Vidxh + v^dyh. 
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(19) 



For normal stress, we assume that the jump across the interface is balanced by surface tension. 
That is, 

k ■ ■ k = 2jH, (20) 

where [cr]" = a — cr^ and cr* = —pol is the stress tensor of the air phase, pq is the atmospheric 
pressure, 7 is the surface tension, H is the mean curvature and k is the surface normal vector 

k = ^ i-d^h, -dyh, if . (21) 

For tangential stress, we assume that there is no jump at the interface 

k ■ ■U = 0. (22) 

That is, we assume that no tangential surface tension gradient exists, as is appropriate for strong 
anchoring. For the case where surface gradient exists, see Ref. 33. 



III. LONG- WAVE HYDRODYNAMIC DESCRIPTION IN TWO DIMENSIONS 

In this section we restrict attention to two space dimensions and focus on the long-wave ap- 
proximation of the governing equations presented previously in Sec. II. The full details are given 
in Appendix A. The aim here is to study the contribution of nematic elasticity to the free surface 
evolution, and to distinguish results obtained using different scalings and boundary conditions. 

Assume the flow is two dimensional and y-independent, so that the director field can be ex- 
pressed as n = (sin^, cos 6')-^ where the angle 9 is taken as the difference between the director 
orientation and the positive z-axis, as shown in Fig. 1(a), and the velocity field is i> = (u, wf. We 
introduce long-wave scalings to nondimensionalize the governing equations. The scalings are 

{x, z) = {Lx,SLz), {u,w) = {Uu,SUw), t = P=^^P, (23) 

U o'^L 

where U is the scale of fluid velocity, 5 = H/L ^ 1 is the ratio between the typical film thickness 
scale, H, and a typical lateral length scale, L. In addition, in order to focus only on the nematic 
elasticity, we approximate the nematic viscous stress tensor by its Newtonian equivalent, setting 
aYj = 2neij, where /x = 0:4/2. 
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A. Weak elasticity 



Assuming that the elastic free energy is weak compared to the pressure, we can introduce the 
dimensionless number (Ericksen number) 

kJJ^. (24, 

The leading order equations are then given by (after dropping the over-bars) 

d^-p = dlu, (25a) 

d,p = 0, (25b) 

Kdle = 0, (25c) 

d^u + d,w = 0. (25d) 

In addition, the leading order boundary conditions are 

9{z = 0) = '^, 9{z = h) = 0, (26) 

P = Pq-C d^h, d,u = 0, atz = h, (27) 

where C = fiU is the inverse capillary number. 

It is easily seen that the velocity field and the director field are decoupled. The film thickness 
evolution equation is obtained as 

dth + dj^h'dlh\ =0, (28) 



and the director field satisfies 



2 V h 

This corresponds to the approach taken in Refs. 25-21. 



= 1(i-t)- (29) 



B. Moderate elasticity 



Instead, if we introduce the Ericksen number as 



/? = ^. (30) 
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the leading order equations are given by (after dropping the over-bars) 




(31a) 



dzP 

Kdle 



0, 



0, 



(31b) 



(31c) 



dxU + dzW 



(31d) 



with leading order boundary conditions 



e{z = 0) 



71 



9{z 



h) = 0, 
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(32) 



p = Po-Cdlh- K{d,ef, -K{djd,9 + {d,efd^h) + d,u = 0, at z = h. (33) 



Under such a scaling, the director field is decoupled from the flow and is given by Eq. (29). 
The tangential stress boundary condition in Eq. (33) is then reduced to dzu{z = h) = 0. We can 
therefore solve for the pressure and velocity field exactly. As a result, the film evolution equation 
is given by 



where K = 7t'^K/4. 

In contrast, Ben Amar & Cummings^ employ Eqs. (31, 32) and impose the normal stress bal- 
ance assuming that the jump of the pressure is balanced by surface tension alone, as is appropriate 
for Newtonian fluids; that is, they use Eq. (27) instead of Eq. (33). As a result, they obtain an 
equation much like Eq. (34) but with the opposite sign for the elasticity term. This issue will be 
discussed later in Sec. V. 

IV. GRADIENT DYNAMICS FORMULATION FOR A THIN FILM OF NEMATIC 
LIQUID CRYSTALS IN TWO DIMENSION 

It was noted some time ago that the time evolution equation for the height of a thin Newtonian 
film on a solid substrate can be written in a variational form in situations where inertia can be 
neglected^ The evolution of the film thickness h follows a dissipative gradient dynamics 
governed by equation 




(34) 




(35) 
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where 5 /5h denotes functional variation with respect to h. The resulting relaxation dynamics is 
governed by the free energy functional F with the mobility function Q{h). 

Such an approach may also be used to obtain the evolution equation for a NLC film in the limit 
of moderate elasticity discussed above in section IIIB. Restricting our attention again to a 2D 
geometry, we simplify the elastic distortion energy for the case of lateral long-wave distortions, 
i.e., we assume the scalings given in Eq. (23). The bulk elastic energy is to leading order 

K 



(36) 



We further assume that the director adjusts instantaneously to its steady state as compared to the 
fluid relaxation time, i.e., we assume K = 0{fiUL) . Then, the director field can be solved for 
exactly to obtain a linear profile as shown in Eq. (29) assuming strong planar anchoring at the 
solid substrate and strong homeotropic anchoring at the free surface. The corresponding director 
orientation across the film is sketched in Fig. 1(b). 

As a result, the bulk elastic energy wp (energy/volume) of the NLC can be rewritten as: 



The free energy functional is then expressed as 



K 
2h?' 



Cds + 



Wp dz ] dx 



C 1 + 



+ hwi 



dx. 



(37) 



(38) 



where ds [1 + {dxhy/2) dx is the approximated surface element. The evolution equation of the 
film is given in gradient dynamics formulation by introducing F into Eq. (35): 

K 



-dr. 



Qih)dr Cdih + 



2h^ 



(39) 



where the mobility function Q{h) can be obtained from the Poiseuille NLC flow, Eq. (A25). One 
should note that Eq. (39) and Eq. (34) are identical when Q{h) = 



A. Linear stability analysis 

To have a basic understanding of the elastic contribution to the stability of NLC free surface, 
we analyse the linear stability of a flat film, h = h^. Assuming h = Hq + ^, ^ <ti in Eq. (39), to 
leading order we have 

dt^ = Qiho)\~Cdt^ + ^^dl^\. (40) 
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With the harmonic mode ansatz ^ = exp{ikx + ujt) one obtains the dispersion relation 

u = -Q{ho)(^k^ + ^^k'^. (41) 

Note that the constants C, K and the film height are always positive and therefore the growth 
rate uj is never positive for any wavenumber k. This implies that the elastic term is always stabilis- 
ing and in the case of strong anchoring the flat film h = ho'is always stable if only capillarity and 
elasticity are taken into account. 



V. DISCUSSION AND CONCLUSION 

We have presented several approaches to derive the evolution equation for free surface films of 
nematic liquid crystals with strong anchoring at both interfaces. We have demonstrated the consis- 
tency between the long-wave approximation model, Eq. (34) and Eq. (A24), and the model derived 
through a thermodynamically motivated gradient dynamics formulation, Eq. (39). Moreover, the 
elastic energy contribution acts in a stabilising manner in each of these models, consistent with the 
physically-motivated arguments of, e.g., Ref. 9 and 36. 

In contrast, the long-wave models of Refs. 7, 22-24, which use an alternative normal stress 
balance, lead to qualitatively different results. The normal stress boundary condition in these 
papers neglects the contribution of the elastic stress tensor, which leads to a change in sign of 
the elastic contribution in the free surface evolution equation. A third approach used by Carou et 
al.^^"^^ scales the nematic elasticity such that, to leading order, the free surface is unaffected by 
the elasticity, and one recovers the Newtonian thin film equation, Eq. (28). 

One should note that the strong anchoring models presented here are only valid for rather thick 
films as noted in Sec. I. First, the main assumption of the model - the strong anchoring of the 
director at both interfaces - is only valid for h ^ he where he is defined in Eq. (1). For moderate 
film thickness h ^ he or even thinner, the surface anchoring energy has to be taken into account. 
This may be done via an ad-hoc amendment of the free surface anchoring condition in Eq. (32) 
(cf. the approach taken in Ref. 22, for the model with the 'Newtonian' normal stress balance). 
Alternatively it may be modelled via the variational approach that will be the subject of future 
work. 

Second, as for isotropic liquids with film thickness below about 100 nm, long- and short-range 
effective intermolecular forces between the substrate and the free surface have to be taken into 
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account possibly through a Derjaguin or disjoining pressure that describes wettability effects . 
For nematic liquid crystals the influence of van der Waals interactions has been discussed, e.g., 
in Ref. 8 and 38. Additional Casimir-type forces may be induced by fluctuations of the director 
orientation, most notably in very thin films with uniform director orientation, i.e., in the planar (P) 
state^^'^^. Note, however, that the notions "disjoining pressure" or "structural disjoining pressure" 
are used in Refs. 40-42 to denote the pressure contribution resulting from the elasticity of the 
liquid crystal, i.e., the last term in Eq. (39). 

We would also like to point out that, within the present long-wave scalings (Eq. (23)), there 
is no distinction in the elastic energy whether the director is bent clockwise or counterclockwise. 
The two director profiles shown in Fig. 1(b) (on the right hand side and on the left hand side of 
the dashed (black) line) have exactly the same elastic energy, K /2h?'. However, such a situation is 
still not allowed even though the elastic energy is continuous across the dashed (black) line. The 
director field is discontinuous and it breaks the long-wave assumption {dl9 <ti dl9). A simple way 
to circumvent this was proposed in Ref. 23, whereby the discontinuity of Fig. 1(b) is smoothed 
out over a given range. More sophisticated models for real defects are needed. For instance, one 
may incorporate a description of the dynamics of the scalar order parameter related to the nematic- 
isotropic transition. Away from the phase transition it can be employed to model defects. Such 
a model would also allow one to tackle the structuring of films that occurs close to the nematic- 
isotropic transition^'' 

In conclusion, we have clarified how the elastic contribution influences the free surface of a 
nematic film under the strong anchoring assumptions. Within the long-wave scalings, we have 
discussed two cases, corresponding to weak and moderate elasticity, respectively: 

• = O (5): The bulk elasticity has only a minor influence on the free surface evolution. 
It does not affect the stability of a film. The evolution of the film and the director field are 
given by Eq. (A18) and Eq. (A20), respectively. 

• = 0(1): The strong antagonistic anchoring makes a significant contribution leading to 
a diffusion-like term in the film surface height evolution equation - see Eq. (39). Further- 
more, the director always maintains its steady state, given by Eq. (29). 

The models can be derived either from asymptotic expansion of the nemato-hydrodynamic equa- 
tions or from a thermodynamically motivated gradient dynamics formulation. It is found that the 
elastic distortion energy is always stabilising. 
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(a)Coordinates (b)Director orientation 

FIG. 1. (Color online) The coordinates used in this manuscript are given in (a). The angle 9 of the director 
is measured with respect to the positive z axis. In (b), we present two possible director profiles of a hybrid 
film. The molecules can bend either clockwise or counterclockwise across the film. The solid (blue) top 
curve indicates the free surface, the solid (black) bottom curve indicates the solid substrate, short (red) lines 
represent the orientation of the director field, and the dashed (black) line indicates the defect location. 
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Appendix A: Long wavelength approximation of a thin film of NLC 

1. Ericksen-Leslie equations in two spatial dimensions 

Assume the flow is two dimensional and y-independent, then the director field can be expressed 
as n = (sin 9, cos 9)^ and the velocity field is v = (u, w)'^. The elastic energy reduces to 

^F = ^ {{djf + {d,9f) . (Al) 
13 



Without fluid inertia, the linear momentum equations are then given by (from Eq. (13b)) 



5r 



-Kdj {die + die) + + d,a 
-Kd,e {die + die) + d^a^^ + d,a^ 



V 
121 



V 
22- 



(A2) 
(A3) 



(The viscous stress tensor, cr^, is defined in Appendix B, Eq. (Bl).) For the angular momentum 
equation, Eq. (13c), one can eliminate the Lagrange multiplier A by performing an inner product 
with the vector = (cos 6*, — sin e)'^. We then have 



K{dle + die) = 71 



+ y [{da^u - d^w) sm{2e) + {d,u + d^w) cos(20)] . 



(A4) 



The continuity equation, Eq. (13a), is rewritten as 



dxU + dzW = 0. 



(A5) 



a. Boundary conditions 



In 2D, the boundary conditions for the director field, assuming strong anchoring at both inter- 
faces (planar at the substrate and homeotropic at the free surface), are 



e{z = 0) 



[z = h) = COS 



-1 



1 



For the velocity field, we assume no-slip at the solid substrate, 

u = w = 0, at z = 0. 
At the free surface we have the kinematic boundary condition 

w = dth + ud^h, at z = h, 
which can be combined with the incompressibility condition, Eq. (A5), to be 

dth + udz ] = 0, 



or equivalently. 



h 

dth + dx [ I dzU {h — z) dz ] = 0. 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 
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(Note that the no-slip boundary condition, u{z = 0) = was imposed in deriving Eq. (AlO).) 

For the balance of normal and tangential stresses, we first note that the stress tensor for a NLC 
film is written as 





1 




cr = —p 






1 





+ 

















(All) 



and the stress tensor of the air phase is cr* = —pol. We assume that the jump in the normal stress 
is balanced by surface tension and the jump in tangential stress is zero. That is 



(A12) 



where k is the curvature, k is the normal vector at the free surface and t is the tangent vector at 
the free surface, defined as 

dlh . 1 . „ , 1 



respectively. 



2^,3/2 • 



k 



-d^h, 1] 



t 



:i,d,hy, (A13) 



2. Non-dimensionalisation and long-wave approximation 

We make the usual long- wave scalings to nondimensionalize the governing equations as shown 
in Eq. (23). Also we rescale the coefficients of nematic viscosity by the Newtonian equivalent, 
setting ai = yua, where /i = 0:4/2. For the elastic constant, we assume K = ejJJ LK where e is a 
parameter of order o{\/5) that will be specified later. 

The leading order equations are then given by (after dropping the over-bars) 

d.^V+^'-^id^ef^ =-eKd,edl9 + d,{qi{9)d,u), (A14a) 

9z{p+'-^ (dzey^ = -eKd^edje, (A14b) 
eKd'^^e = -6q2{e) d,u, (A14c) 
d^u + d,w = 0, (A14d) 

where qi{9) and q2iO) are related to the viscous stress tensor, their full expressions are given later 
in Sec. B. The leading order boundary conditions are the kinematic boundary condition, Eq. (AlO), 
with 

e{z = 0) = -, e{z = h) = 0, (A15) 
2 
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p = p^-C dlh-eK{d,ef, -eK{dJ d,9 + {d,9f d^h) +qi{9) d,u = 0, at z = h, (A16) 
where C = / jjiU is the inverse capillary number. 



a. Weak elasticity (e = 6) 

Assuming the elastic free energy is weak compared to the pressure, we can choose e = 5. 
Observing that Eq. (A 14b) reduces to = at leading order, one can solve the pressure exactly 
and the velocity is then determined by 

d,u{x, z) = {h - z) dlh. (A17) 

Hence, by using Eq. (A 10), we obtain the film evolution equation as 

ht + Cd^{Q{h)dlh) (A18) 

where 

Q{h) = [ ^^^^ dz. (A19) 

Jo qi{9) 

In addition, the director field satisfies 

dle = ^''^^(z-h)dlh (A20) 

with boundary conditions defined in Eq. (A 15). 

One can see that the nematic elasticity as well as viscosity only have influence on the mobility 
function Q, and thus have no influence on the stability of a free surface. This formulation has been 
studied extensively by Carou et alP'^^ both analytically and numerically under the assumption of 
small director variation. 

b. Moderate elasticity (t = l) 

On the other hand, if we have e = 1, Eq. (A14c) reduces to dl6 = at leading order and hence 
the director field reaches a linear profile in z as shown in Eq. (29). Moreover, Eqs. (A 14a) and 
(A 14b) are simplified to 

5. I p + j = 9, (gi(0) d^u) , (A21a) 
d,p = 0, (A21b) 
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with boundary conditions at the free surface 



K 

p = Po-C dlh - — , = 0. 



We can therefore solve the pressure and velocity field as 



K 

p{x,z) =po-Cdlh - —, d;,u{x,z) = 
As a result, the film evolution equation is given by 



K 



ht + 



Qih) lCdih--d,h 



-Cdth+—d.,h 



K 



where Q{h) can be evaluated explicitly as' 



T^J Jo QliO 



z — h 



{All) 



(A23) 



(A24) 



(A25) 



Appendix B: Viscous stress tensor of nematic liquid crystal in two dimension 



The viscous stress tensor of NLC in 2D is written as 



.V 



ai (sin^ 6 d^u + cos^ 6 d^w + sin 6 cos 6{dzU + dxw)^ 



sin^ 6 sin 6 cos 6 
sin 9 cos 9 cos^ 9 



' Di^ 2 



+"3 



sin 9 cos 9 cos^ 9 

— sin^ 9 — sin 9 cos 9 

sin 9 cos 9 — sin^ 9 

— cos^ 9 — sin 9 cos 9 



+ 



+ 



+ 



y 

2 

"6 



(Bl) 



Dt 2 

2c}2;M dzU + (9a;ty 

(9j;M + dxW 2dzW 



2 sin^ 6' dxU + sin cos 9{dzU + S^^w) 2 sin 9 cos 6' S^jM + cos^ 9{dzU + 

sin^ 9{dzU + ^^.w) + 2 sin 6* cos 9 dzW sin cos 9{dzU + ^^.w) + 2 cos^ 9 dzW 

2 sin^ 9 dxU + sin 6^ cos 9{dzU + sin^ 9{dzU + S^^w) + 2 sin 6' cos 9 dzW 

cos^ 6'(92'U + dxw) + 2 sin cos 9 dxU sin 9 cos ^(S^m + dxW) + 2 cos^ 6^ dzW 



Similarly, the coupling term, g, between the director and velocity field can be written as 



-7i 



dzU — dxW 
2 



cos 9 
— sin 6 



72 

2 



2 sin 6^ dxU + (S^m + ^^..w) cos 9 
{dzU + dxw) sin 6* + 2 cos 9 dzW 



(Bl) 
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We also note that, within the long- wave scalings [Eq. (23)], to leading order we have 

ai; = mm d.u + O (^^^ , n^ g = -^2(9) d^u + O (^^^ , (B3) 

where n = 04/2, n-*- = (cos6', — sin^^)^ and 

Wi(^) = 2 ["^4 + ^'^1 ^^^^ ^ '^'-'^^ ^ ('^5 ~ '^2) cos^ 6 + (as + as) sin^ 6*] , (B4) 
Hq2{e) = ^[li-l2COs{2e)]. (B5) 

As an example, for Newtonian fluids, qi{6) = 1 and ^2(6*) = 0. 



Appendix C: Notation conventions 

For clarity, we list all the notations used. We write for a vector n = rii ox m = nii, for a 
tensor cr = aij or k, = Kij; as the superscript T denotes transposition one has cr^ = aji. Further, 
eijk is the alternator. The notations for operators and products are Vn = j, V ■ n = fc, 
V X n = eukUk^i, V cr = a,k,k, cr ■ k = aikK^j, cr : k, = dkiKik, cr • n = aikUk, nxm = eukUkmi, 
where ', z' denotes the partial derivative with respect to the ith component. 
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